|^lb-AOM  555 
UNCLASSIFIED 


NAVAL  UNDERWATER  SYSTEMS  CENTER  NEW  LONDON  CT  NEW  LO— ETC  F/G  17/1 
A  NEARFIELD  MODEL  OF  THE  PARAMETRIC  RADIATOR  PART  II.  SATURATED— ETC  CUI 
NAY  76  R  H  MELLEN 
NUSC-TN-PAA-53-76 


NL 


AO  A  i)  B  3  5  5  5 


t 


% 

A  Y«¥ 


i 


TM  No. 
^EAU-53-76 


i/ 


NAVAL  UNDERWATER  SYSTEMS  CENTER 
HEW  LOUDON  IABORATOHT 
HEW  LONDON,  CONNECTICUT  06320 


/  V  \  Technical  Memorandum 

v  /  J  ^  i 

),  A  NEAHFTEID  MpDEL  OP  THE 

JART  Ha  SATURATED  1SJURCE 


IATOR 


•  /  '  '*17  Msy  *76  ' 


Prepared  Tax  ^  ^  — _ 

/  H.^feHen 

’'  Hans  and  Analysis  Directorate 


!  j L:  I  l\ }  y  O  -• 


?  ! 


c 


tffeELEC 

W  apo  ~  p, 

r:: 

; 

',930  :  '  ■ 

u 

A 

■1-w 

( 


♦To  be  presented  at  the  7th  International  Symposium  on  Nonlinear 
Acoustics,  Virginia  Polytechnic  University,  Blacksburg,  VA, 
August  1976* 


kj. 


****** 


cl4(t 


Reverse  Blank 

Approved  for  public  release;distribution  unlimited. 

80  4  28  009 


/  ■"  'J 


HOCEDM0  PAGE  BUNK-NOT  HIKED 


TM  No. 
PA4-53-76 


ABSTRACT 


vl/ 


In  the  ne airfield  model  of  the  parametric  radiator,  source 
levels  were  obtained  by  assuming  the  primary  wave  to  be  plane- 
cylindrical  or  diverging-conical ,  depending  on  whether  the 
observation  range  was  less  or  greater  than  the  Rayleigh  distance. 
The  two  cases  are  now  combined  by  means  of  the  "horn "  model, 
which  approximates  primary  wave  diffraction  for  a  piston  source. 
Saturation  effects  are  also  incorporated  by  means  of  a  nonlinear 
attenuation  approximation. 
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INTRODUCTION 

la  reference  (1),  a  coupler  contour  integration  method  was  used 
to  calculate  nearfield  parametric  source  levels  for  two  primary  wave 
models:  a  plane-wave  cyllndrically  collimated  source  for  ranges 
less  than  the  Rayleigh  distance  and  a  spherical  wave  conically  colli¬ 
mated  source  for  greater  ranges.  The  two  cases  are  now  combined  by 
means  of  the  horn  model  and  the  effects  of  saturation  are  included 
by  means  of  a  nonlinear  taper  approximation. 


BOHR  MODEL 


m  the  horn  model  of  Fig.  1,  we  require  that  the  primary  wave 
be  predominantly  plane  near  the  source  and  predominantly  spherical 
for  ranges  much  greater  than  the  Rayleigh  distance  Rq  .  The  cross- 
sectional  area  of  the  horn  is  chosen  to  maintain  constant  source 
density  per-unit -length  when  attenuation  is  neglected.  On  the  axis, 
the  primary  peak  pressure  is  approximated  as 


2 

p  — . 

P  -  Re  YToZ  ^  «cp  (ia*t) 
i  n»l 


(1) 


where  PQ  is  the  initial  pressure  of  each  component,  ?  is  the 
retarded  time  and  f  ■  r/Ro  is  the  scaled  distance  from  the  piston 
source.  J 

Squaring  Eq.  (l)  and  taking  the  cross-sectional  area  as 
S  -  Sq  (l  +  /■  2)  «e  can  write  the  dimensionless  normalized  source 
density  per-unit -length  as 
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2  (-CVQ  +  f  2  cos2  (ilr)] 

1  2 


m  1 


(2) 


where  Sq  is  the  piston  area,  -CL  ■  +  '**2)/2  and  the  bars 

indicate  the  average  value.  The  sin2  (At)  term  is  taken  as  the 
cylindrical  component  of  the  source  density  and  the  cos2  (117) 
term  as  the  conical  component. 
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The  parametric  source  strength  relative  to  the  We  starve  It  source 
strength  at  infinite  range  can  then  be  written  symbolically  as 


T  dTv  +  T  r2*1 

<>  ITT?  '-'o  i  +f 


where  Yi  and  Y2  are  Eq.  (9)  and  Eq.  (12)  of  reference  (l)  * 
respectively. 


SATURATION  TAPER 

Reference  (2)  gives  the  derivation  of  a  nonlinear  taper  function 


T2(r) 


K*L  /I  +  2  K 


where  <T  (r)  is  the  fictitious  nonlinear  phase  distortion  angle  at 
the  peak  envelope  of  the  primary  wave.  For  <T  >  1  ,  the  waveform 
becomes  multivalued  and  the  jump  approximation  is  used  to  restore 
continuity  at  the  shock  front.  Interaction  effects  of  absorption 
are  ignored  and  the  linear  and  nonlinear  attenuations  are  multi¬ 
plicative. 

The  phase  distortion  angle  is  taken  as  the  magnitude 


<T(r).|ta  /  /F<r-)l  dr' 

jo  Co2  ^  o 


where  £  is  the  nonlinearity  number  and  |F(r')l  is  the  peak  envelope 
pressure  of  Eq.  (l).  Equation  (5)  becomes 


<T  (j*  )  «  2  X 


/  _  a  f' 

./  1  j.  JT.* 


-  2  X  sink"1  (P  ) 


TM  Bo. 
PA4-53-76 


where  X  «  ft  P0  ko  W/o  co2  ia  the  saturation  number. 
Equation  (3)  then  becomes 


vJ  P  +  r  «2/2 £(f) 

*4  “77p“  Jo  i  +  f  2 


(7) 
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For  the  cylindrical  component  Eq.  (9)  of  reference  (1)  can  be 
written 


1 

i  TL 


J°  dt  a**t  -  |  exp  (-  L  v^)  J~ 


oo 


(8) 


where  L  is  the  range  from  the  origin  to  the  field  point,  Uo  -  2<*L  , 
t-2str,  to  -  2* Ho  ,  n-*0k/koUo,  B1  *  V/*0  T1  1  “J 

z.  «  t  -Uq  +  ^  TL  •  The  oecoixl  integrand  of  Eq.  (8j  involved  a  change 
to  the  complex  variable  z  ■  2*  *  *  *  k  (  ^a*  +  x^  +  x)  where 
x  •  r-L  and  a  change  of  contour  to  one  of  constant  imaginary  value. 

In  Eq.  (6)  and  Eq.  (7),  the  corresponding  variable  f  becomes  complex 

ew3  hovA 


and  we  have 


<1  ^ 
1"  L  to  ko 


Uo 
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dt  e~^  T2  (t/tp)  i  e 

1  +  (t/to)k  ? 


oe  .  jL-  1 

-  ‘  Tl  J*  dt  e-t (l-  v/*l2+Bl5)  ^(Jl)  ( 
o  lT??  Ii 


(9) 


where  ^fl  »  +  2u0-  i/2?  +  ®12  )/2  to 
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For  the  conical  component,  Eq.  (12)  of  reference  (l)  can  be 
written 


Uq  w 

Vo  J*  dt  e^  J"  dt  e-t  (1-  /zg2  +  ®22) 


+  2uq  “  /z2z  +  B22 


where  t2  -  Uq  k/ko  to  ,  B22  *  **■  *  Uq  V2  and  z2  »  t  -Uq  +  i  v2  . 

The  residue  is  omitted  because  the  combined  solution  has  no  pole. 

By  «•<"«<  change  of  ^  to  complex  form  in  Eq,.  (6)  and  Eq.  (7), 
we  have 

Uo  00  *2 

%  -  f  dt  e-*  T^t/tp)  if  dt  e-t  J2  (1.  yzg^Bg2 )  T2^) 

c  150  O  1  +  (t/t0)2  ?  O  1  +  ^22 

_  (11) 

when  52  *(z2  +  2Uq  -  /z22  +  ?22^/2  *o  » 

Equations  (9)  and  (11)  were  evaluated  numerically  as  discontinuous 
functions  at  t  ■  Uq  . 

HESUIiTS 

Some  results  are  shown  in  Figs.  2-17  for  selected  values  of  the 
downshift  ratio  k^k  ■  F0/F  ,  the  absorption  number  AFo(dB)  * 

8.7oS>Ro  ,  and  the  saturation  parameter  X(dB)  *  20  logio  X  .  The 
ordinates  are  Y(dB)  »  20  log  1 Y  |  and  the  abscissas  are  the  scaled 
range  L/?q  .  For  values  of  the  saturation  parameter  less  than  -1 5 
dB,  the  effects  of  saturation  rapidly  become  negligible.  At  long 
ranges  the  relative  levels  approach  -20  log™  X  dependence  when  X 
is  large.  Since  the  reference  source  level  (Westervelt )  is  given 
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by  20  login  *o  Y0  ■  PoBoX(F/F0)2/4<CRo  >  the  actual  source 

levels  go  from  Uo  logxo  X  dependence  in  the  unsaturated  regime 
to  20  logxo  X  dependence  in  the  saturated. 

The  RMS  parametric  source  level  nay  be  calculated  from  the 
formula 

SL  -  SLq  +  Y(dB)  +  X(cLB)  -  40  log^  (F0/f)  -  20  log  (ARq) 

(12) 

+  6.7  dB//l|iPa  •  m 


where  SLq  is  the  RMS  source  level  of  each  primary  component.  For 
10 °C  sea  water  at  moderate  depth  we  have  from  reference  (3) 


X(dB)  -  31^  +  20  logio  (FqA  Mz)  -  280.6  (13) 


where  F0  is  the  mean  primary  frequency  in  kHz. 
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RELATIVE  SOURCE  LEVEL  VS  SCALED  RANGE 
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RELATIVE  SOURCE  LEVEL  VS  SCALED  RANGE 
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RELATIVE  SOURCE  LEVEL  VS  SCALED  RANGE 
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